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of Planck scale suppressed Lorentz Invariance Violation (LIV) on the propagation of 
cosmic ray nuclei. In particular we focus on LIV effects on the photodisintegration of 
nuclei onto the background radiation fields. After a general discussion of the behavior 
of the relevant quantities, we apply our formalism to a simplified model where the LIV 
parameters of the various nuclei are assumed to kinematically result from a single LIV 
parameter for the constituent nucleons, t], and we derive constraints on rj. Assuming 
a nucleus of a particular species to be actually present at 10^° eV the following 
constraints can be placed: —3 x 10~^ < 77 < 4 for ^^Fe, — 2 x 10"^ ^ ^7 ^ 3 x 10"^ 
for ^^O and —7 x 10~^ ^ ''7 ^ 1 x 10^^ for ^He, respectively. 
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1. Introduction 

Possible small deviations from the exact local Lorentz Invariance (LI) of general 
relativity have received a growing interest in recent years. On the theoretical side, 
hints of Lorentz Invariance Violation (LIV) were found in various approaches to 
Quantum Gravity (QG) [1, 2, 3]. On the observational side, high energy astrophysics 
observations allowed to probe the consequences of LIV in various contexts [4, 3]. 

Lorentz symmetry breaking is not a necessary feature of QG, but it is clear that 
any possible LIV effect connected with the Planck scale could provide an observa- 
tional window into QG. However, to directly observe phenomena connected with the 
Planck scale Mpi = 1.22 x 10^® eV would require the center of mass energy of, e.g., 
a scattering process to be comparable to Mpi. This is 15 orders of magnitude larger 
than what the LHC can probe with its design center of mass energy of 14 TeV. On 
the other hand, if we are testing LI specifically, then also non-LI quantities can be 
important. The energy of the particle in some frame, or a cosmological propaga- 
tion distance are widely discussed examples. These quantities can be so large as to 
effectively offset the Mpi suppression to a physical observable, so that very small 
corrections are magnified. For this reason, they are called "windows on QG" . 
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In order to correctly identify such "windows on QG" it is important to place them 
into a dynamical framework. A standard method is to study, within the context of 
EFT, a Lagrangian containing the standard model fields and all LIV operators of 
interest that can be constructed by coupling the standard model fields to new LIV 
tensor fields that have non-zero vacuum expectation values [5, 6, 7, 8, 9, 10]. 

A word of caution should be made at this point. Lorentz invariance of physical 
laws relies on only few assumptions: the principle of relativity, stating the equiva- 
lence of physical laws for non-accelerated observers, isotropy (no preferred direction) 
and homogeneity (no preferred location) of space-time, and a notion of precausality, 
requiring that the time ordering of co-local events in one reference frame be pre- 
served [11, 12, 13]. The above described procedure leads to breaking of the principle 
of relativity, but one may wonder if other possibilities exist that violate Lorentz in- 
variance but preserve the principle of relativity. One such possibility is represented 
by the very special relativity framework [14], which corresponds to the break down of 
isotropy and is described by a Finslerian-type geometry [15, 16, 17]. In this example, 
however, the new relativity group generators number less than the ten generators 
associated with Poincare invariance. An example of a new relativity group different 
from Poincare but with exactly ten generators is not known in commutative space- 
times. In non-commutative spacetimes, however, such a construction is possible and 
was termed "doubly" or "deformed" (to stress the fact that it still has 10 genera- 
tors) special relativity, DSR [18]. Unfortunately, the various DSR candidates face in 
general major problems regarding their physical interpretation (they do not yield a 
low energy EFT, in general) and a working model is not yet available (see however 
[18] for recent attempts in new directions). 

We now turn back to our preferred EFT approach. A generic result of the pro- 
cedure described above is the presence of modified dispersion relations for particles, 
of the form 



where m is the particle mass, E its energy, the function / represents the QG con- 
tribution and can depend generically on the momentum p, on Mpi and on some 
intermediate mass scale fi. For simplicity we assume that only boost invariance is 
broken, while rotations are preserved (see [4] for further comments on rotation break- 
ing), so that / depends on p = \p\, rather than on p. Moreover, at <^ Mpi we can 
expand / so that eq. (1.1) reads 



where r/^"^ are constant dimensionless coefficients to be constrained. Terms with 
n = lead to strong deviations at small energies, while n > corresponds to high 



E^-p^ = m^ + fip,Mpufx), 



(1.1) 




(1.2) 



-2- 



energy deviations from Lorentz symmetry. The latter situation is what we will discuss 
in this work. 

Many of the parametrized LIV operators have been very tightly constrained via 
direct observations (see [4, 1, 2, 3] for reviews). In particular, terms n = 2, coming 
from dimension five and six CPT even LIV operators, have been recently directly^ 
constrained [20] in the hadronic sector by exploiting ultra high energy cosmic ray 
observations performed by the Pierre Auger Observatory (PAO) [21]. Indeed, the 
successful operation of the PAO has brought UHECRs to the interest of a wide com- 
munity of scientists and already allowed to test fundamental physics (in particular 
Lorentz invariance in the QED sector) with unprecedented precision [22, 23, 24]. 

The UHECR constraints [25, 26, 27, 28, 29, 30, 31, 20, 32] rely on the behavior 
of particle reaction thresholds with LIV. What matters for threshold reactions in the 
presence of modified dispersion relations as in eq. (1.2) is not the size of the LIV 
correction compared to the absolute energy of the particle, but rather the size of the 
LIV correction to the mass of the particles in the reaction. Hence the LIV terms 
usually become important when their size becomes comparable to the mass of the 
heaviest particle. This criterion sets the presence of a critical energy E^r above which 
LIV effects are relevant in a given threshold reaction. If the LIV term scales with 
energy as E"'~^^, then Ecr ~ (m^Mpj)"+2 [33]. According to this reasoning, the larger 
the particle mass the higher is the energy at which threshold LIV effects come into 
play. 

UHECR constraints have relied so far on the hypothesis, not in contrast with any 
previous experimental evidence, that protons constituted the majority of UHECRs 
above 10^^ eV. Recent PAO observations [34], however, showed strong hints of an 
increase of the average mass composition with rising energies up to E' ~ 10^^'^ eV, 
although still with large uncertainties. Hence, experimental data suggests that heavy 
nuclei can possibly account for a substantial fraction of UHECR on Earth. Previous 
constraints need then to be revised. The purpose of the present work is to study 
the effects of LIV in the propagation of UHECR nuclei to establish more firm, albeit 
perhaps less stringent constraints. 

This paper is composed as follows: in sec. 2 we discuss our LIV framework, while 
in sec. 3 we present general facts about propagation of UHECRs in the intergalactic 
medium. In sec. 4 we deal with photodisintegration, the main process affecting the 
propagation of nuclear species, and develop a formalism to take possible LIV effects 
into account. The results of these calculations and some of their implications are 
finally presented in sees. 5 and 6. 

^Notice that all these operators can be indirectly constrained by EFT arguments [19], as higher 
dimension LIV operators induce large renormalizable ones if we assume no other relevant physics 
enters between the TeV and Mpi energies. Supersymmetry, however, is an example of new relevant 
physics that can change this argument. 
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2. Theoretical Framework 



In order to study the phenomenological consequences of LIV induced by QG, the 
existence of a dynamical framework in which to compute reaction thresholds and 
rates is essential. We assume that the low energy effects of LIV can be parametrized 
in terms of a local EFT and, for simplicity, that rotation invariance be preserved. 
Therefore we introduce LIV by coupling standard model fields to a non-zero time-like 
vector. 

We take the opportunity offered by the high energy involved in UHECR physics 
to focus our discussion on terms suppressed by two powers of Mpi. The general 
discussion of this case can be found in [9, 10], where the structure of the LIV terms 
and the full set of modified dispersion relations can be found. Accordingly, the final 
dispersion relation for nuclei we want to further elaborate on is 

E'=p' + m' + v('^^, (2.1) 

with rj independent of the helicity state of the nucleus, having assumed parity to 
be preserved. This is just a simplifying assumption that reduces the number of 
degrees of freedom we have to deal with in the following. Having already extracted 
the relevant powers of Mpi, r] should be of order 0(1) if the LIV effects are due to 
physics at the Planck scale. Therefore, any limit |?7| ^ 1 is effectively a constraint 
on the assumption that QG physics occurs at the Planck scale. 

2.1 Modified Dispersion Relations for heavy nuclei 

The modified dispersion relations (MDRs) derived according to [9, 10] are valid for 
elementary particles, such as electrons, neutrinos and quarks. Other MDRs hold for 
gluons, and the final MDR of a compound particle like a proton is a combination 
of all the LIV of the constituents. According to this approach [35], a parton model 
needs to be assumed for protons and the net proton LIV is determined by the LIV 
terms for the partons along with the parton fraction at UHECR energies. However, 
we will not establish constraints on the bare parameters here, nor are we interested 
specifically in proton dispersion relations. Rather, we assume that each individual 
nucleus has its own independent MDR. 

Therefore, a nucleus with mass A and charge Z will have the following MDR 

Elz = pXz + mlz + Vf/^ , (2.2) 

where Ea,z is the energy, pA,z is the absolute value of the 3-momentum, mA,z the 
mass and rj\ ^ the Lorentz- violating parameter of the nucleus. We assume that only 
nuclei have MDRs. Indeed, the target photons have too low energy for the LIV 
effects to be relevant in their dispersion relations. 
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Now we make a further simplification. Assuming a different rj^^z for each differ- 
ent nuclear species and isotope would lead to O(IOO) free parameters to constrain, 
making it extremely hard to draw meaningful conclusions from the analysis of photo- 
disintegration patterns. However, we can achieve a huge simplification by assuming 
that energy and momentum of the nucleus are the sum of energies and momenta of 
its constituents [33]. Using the fact that a nucleus is made of protons and neutrons 
and that protons and neutrons have almost the same mass, mi, we have that the 
total energy and momentum of a nucleus {A, Z) are given by Ea = Yli ~ AEi 
and Pa = ^iPi ~ ^Pii respectively, where Ei is the energy and pi is the momentum 
of the particular nucleon. With this approximation, the dispersion relation can be 
written as 



El = {AE,f = {Ap,f + {Am,f + = vlz + + • (2-3) 



So now we have only one free parameter, 77, for a nucleon, while for nuclei there are 
effective parameters of the form rjA = rj/A^. Although this is just a phenomenological 
model, we remark that it guarantees that the correct dispersion relations are recov- 
ered when dealing with macroscopic objects [33], where the effects of QG physics 
should be suppressed. It might seem that the problem is now oversimplified. We 
remark, however, that the phenomenology of this model is already very rich for a 
first study of LIV effects on UHECR nuclei propagation, as it will be discussed later. 
We will instead leave to future work the study of a more general parameter space. 

3. Ultra High Energy Cosmic Rays 

The spectrum of Cosmic Rays (CRs) spans more than 12 decades in energy (ranging 
from < 100 MeV to > 10^° eV) with an impressively regular power-law shape 



where the exponent p{E) can be taken as a piecewise function, depending very little 
on the energy [36]. 

One of the most puzzling problems of CR physics is related to their origin. Being 
charged particles, they are deviated while propagating in intergalactic and galactic 
magnetic fields, so that in general their arrival direction on Earth does not point back 
to their sources. As a consequence, the observed CR arrival directions are almost 
isotropically distributed. Only CRs with sufficiently high energy, E > 10^° eV, are 
not deflected significantly in the nC intergalactic magnetic field, leading to a possible 
anisotropy, which was actually detected [37, 38]. 



PI 



dN 
dE 
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The Larmor radius of a relativistic particle with energy E, charge Z and propa- 
gating in a magnetic field B is 

E 1.08 f E \ f B / X 

Accordingly, the fiG galactic magnetic field is able to confine particles within the 
galactic disk up to ~ Z x 10^^ eV, while CRs with higher energies cannot be 
confined and would produce a large anisotropy if they were accelerated in galactic 
source. Therefore, they are thought to have extragalactic origin. 

The maximal energy of galactic CRs being dependent on the charge, we expect 
the actual end-point to the galactic CR spectrum to be determined by heavy nuclei, 
with the chemical composition becoming heavier and heavier with increasing energy. 
Experimental data in the energy range 10^'^ eV ^ E < 10^® eV are compatible with 
this expectation [39]. At energy E > 10^*^ eV, however, the chemical composition 
becomes again very light, being compatible with the presence of just protons [40, 41]. 

Recently, the PAO measured the UHECR chemical composition up to ~ 10^^'^ eV 
[34]. In apparent contrast with findings from other experiments [40], PAO results 
hint at the possibility that a significant fraction of UHECRs with energy E > 10^^ eV 
might be constituted of heavy nuclei. Further observations which may confirm these 
results have been made by the Yakutsk EAS Array [42] . Because of limited statistics, 
it is not clear yet whether nuclei dominate the UHECR fiux also above 10^^'^ eV. 
In the light of these results, it is then mandatory to study the propagation of heavy 
nuclei in the intergalactic medium (IGM). The IGM gas density being very small, 
hadronic interactions in the IGM can be neglected and we can focus on interactions 
of heavy nuclei with photon background fields, in particular the CMB and infrared 
fields. In this respect, the main interaction channels are the following [43]: 

Compton interactions Inverse Compton scattering can be described as the pro- 
cess ^-^+ 7 ~^ 7') the UHE nucleus scatters on a (low energy) photon 
and transfers some of its energy to it. The energy loss is quite small but it is 
remarkable that this reaction, in comparison to the others, has no threshold. 

Electron-positron pair production Pair production, ^iV-l- 7 — )■ ^A^+ e~ -|- e~^, 
only dominates for nucleus energy below 10^° eV [43]. In the case of black body 
radiation at temperature T the inverse attenuation length can be written as 

, _ ^cLE _ ar^Z^ {mekBTf f me \ , . 

Here, a is the fine-structure constant, T is the temperature of the corresponding 
photon field, rg is the classical radius of an electron and rrie its mass. The 
function / is given in [44]. Since the reaction occurs frequently throughout the 
propagation and the energy loss in one interaction is of the order me/nip ^ 
10~^, it may be treated as a continuous energy loss process. 
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Pion production Pion production is a pliotohadronic process of tlie basic type 
iN + -1 N' + X + 7r°'±. Its cross section for single nucleons is well-known 
and can be found in [45]. For nuclei usually a simple model for a combined 
cross section, 

(Ta,z = Zap + (A - Z)cr„, (3.4) 

(Tp and (T„ being the cross sections for a proton and a neutron, respectively, 
is used. In the LI case it begins to become important in comparison to other 
energy losses for energies which lie beyond the energy scales discussed here. 
In some sense, this process can be viewed as a kind of photodisintegration, in 
which the energy available is large enough as to trigger the formation of pions. 

Photodisintegration This is the most important interaction channel for changes 
of the chemical composition, also for the LIV case, hence we will describe it in 
more detail in sec. 4. 

Photodisintegration is the most relevant process for the propagation of heavy 
nuclei above 10^^ eV [43]. In principle, also LIV effects on Compton scattering and 
on pair production should be considered. However, given that these processes are 
only relevant at low energy, we do not expect strong constraints coming from their 
analysis. Therefore, we will neglect them and we will focus on photodisintegration 
in the following. 

The radiation fields relevant for our analysis are the Cosmic Microwave and 
Infrared Backgrounds (CMB and CIB). The ambient photon density per photon 
energy is given by 

dVde vr^ ^kpr _ ^ ' 

T = 2.725 K, for the CMB, i.e. it follows Planck's law of black-body radiation for the 
corresponding temperature. In the case of Cosmic Infrared Background the energy 
dependence is more complicated [46]. 



4. Photodisintegration of UHE Nuclei and its LIV Modifica- 
tions 

Photodisintegration is the main process causing a change of a nucleus species during 
the propagation in the intergalactic medium. It may be described as a process con- 
sisting of two steps: Photo-absorption by the nucleus forming an excited compound 
state which then decays with the emission of (at least) one nucleon [47]. 

The process occurs with the so called Giant Dipole Resonance (GDR), for photon 
energies up to around 30 MeV in the rest frame of the nucleus, while at higher energies 
the quasi-deuteron process, baryon resonances and finally, at extremely high energies. 
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(multi)fragmentation, i.e. the decay of the initial nucleus into several large fragments, 
dominates [48]. 

We are now going to consider a reaction of the form 

iN + ^ ^ ^^,N' + (4.1) 
where A' = A - B and Z' = Z -W . 



4.1 Threshold Computations 

One of the most well known consequences of LIV in MDRs is to modify the threshold 
structure of scattering processes. 

For the purpose of this threshold computation, we will consider, instead of an 
MDR hke in eq. (2.3), a generic MDR hke 

n+2 

7?2 _ ^2 I ™2 I in) ^A.Z r,\ 

^A,z -PA,z + ^A,z + VA,z^rj:;^ ^ (4.2) 

where the case of eq. (2.3) is recovered by setting n = 2 and r]^^\ = rj/A^. 
The energy-momentum conservation equation in the threshold configuration [33] 
reads 



{Ea,z + ey - {pA,z + P7) = iEA',z' + Eb,w) - {PA',z' + Pb,w) , (4.3) 

where Eji,z and the other energy values are given by eq. (2.2) while e and are 
the energy and the momentum of the background photon. The threshold conditions 
are head-on collision of the incoming particles and parallel outgoing momenta [49], 
i.e. PA,zP-y = -PA,zP-f = -PA,ze eiiadpA',z'PB,w = Pa',z'Pb,w- Due to the large energy 
of the nuclei the products are emitted almost parallel to the direction of the initial 
nucleus, i.e. momentum conservation requires p = pa,z ~ Pa',z'+Pb,w ~ {^~y)p+yp 
(0 < y < 1), with y being the inelasticity of the reaction. After some algebra, the 
main threshold equation reads 

("•) _ W /I _ \n+l _ (n) n+1 / 2 2 

'Ia,Z 'IA',Z'\^ y) 'lB,wy n+2 a I^A',Z' ^B,W 



(4.4) 

To find a lower threshold, eq. (4.4) has to be solved for p and minimized with respect 
to y. However, in the case of LIV eq. (4.4) may have two real and positive solutions, 
in which case there are two thresholds - a lower (p^^^) and an upper (p^/i^) one (which 
has to be found by maximizing p with respect to y) [49]. If instead there are no real 
and positive solutions, then the reaction cannot take place at all. Dividing eq. (4.4) 



-8 - 



G 




Figure 1: G{x,y) for n = 2 and different values of 02- The points of intersection with 
the X axis are the thresholds. The following configurations are given: For a2 = ^ and 
02 = (LI) there is one (lower) threshold; a2 = — ^ gives a lower and an upper threshold; 
a2 = — ^ is the tangential value 0^2"^^ ^ ^2 — < Oj""^ the reaction is kinematically 
forbidden. 



by its last term and substituting 

X = — y P, 

"^A',Z' , "^l,W _ ^2 
1-y + y ^A,Z 

we obtain a polynomial of the order n + 2: 

G{x,y) ■.= aniy)x''^'^ + x -1=0, (4.6) 



shown in fig. 1. The number of solutions now depends on a„(?/). 

The region in which the process is generically allowed to take place in some energy 
window is limited by the tangential configuration (G(x, y) = and d^G^x, y) = for 
X > 0), given by 

tang ^ (^ + 1)''+^ .^^g ^ n + 2 

(n + 2)-+2' n + 1 ^ ' 

Using this and the structure of a„(y), we obtain the results summarized in Tab. 1. 
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^A.z ^ Vb,w 
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Ta 
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z ~ Vb,w 
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(n) (n) 

Va,z ^ Vb,w 


in) (n) 

Va,z ^ Va',z' 


PD constr. 


PD constr. 


PD constr. 
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PD constr. 


PD constr. 


excl. 


excl. 




(n) (rt) PI 

Va,z — Va',z' ~ ^ 


PD constr. 


excl. 


LI 


excl. 




> 




excl. 


excl. 






(n) (n) 













Table 1: For some combination of parameters from the region labeled by "PD constr." 
a kinematical configuration for photodisintegration to take place cannot be found so the 
reaction is forbidden and therefore it is possible to set constraints. Inside the white region 
the reaction is always possible at some energy, while the cases marked by "excl." are 

(n) 

excluded since the corresponding combination of conditions for ry^ ^ is not possible. 



In general, constraints on the parameters can be placed by imposing that the 
reaction be allowed and happen between some lower and upper threshold energy 
which are compatible with experimental results. 

4.2 Propagation Lengths 

There are two main quantities to characterize the effects of photodisintegration on 
the propagation of UHE nuclei. On the one hand we have the mean free path (the 
distance between two successive reactions) [50] 

(f) = ^ i.,.^ ? J, * ' ^.'-■-wJ ' • 

with e' = — - — (1 — B)e H V^a\^ and e = e + 28 — - — e. On the other hand 

we consider the attenuation length, i.e. the mean distance over which the particle 
energy is reduced by a factor 1/e with respect to its initial value, ^ 

^ 1 dE dedN{e) f'J, 1 h , . , . 

(4.9) 

In equations (4.8) and (4.9) ethr is the threshold energy for the background photon 
obtained from eq. (4.4), dN{e)/dVde is the ambient photon density, y{e',p) is the 

^It is not possible to define properly an attenuation length for a process that implies the dis- 
ruption of the initial particle, such as photodisintegration is. Here we assume that we can group 
different nuclear species in such a way that the particle "identity" is not changed after a few Xmfp- 
This is reasonable given the poor experimental accuracy in the determination of the mass number 
and charge of an UHECR nucleus. Moreover, assuming that only one or two nucleons are removed 
and that the nucleus Lorentz factor is almost constant, the average cross section changes smoothly 
with the mass and the continuous mass loss approximation can be valid [48]. 
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inelasticity as computed following [51] and f3 = dE/dp. The cross section cr(e') 
used here is taken from [43], but modified by shifting the cutoff to e'^„^ = 10 GeV. 
The integration limits in equations (4.8) and (4.9) are computed following the same 
procedure as [20, 51]. 

For the single parameter model described by eq. (2.3), y{e' ,p) is almost inde- 
pendent of rj. Therefore it is enough to consider only the mean free path Xmfp of 
eq. (4.8) from which the corresponding attenuation length can be obtained dividing 
by ythr,Li = TTT'B,w/ {^A',z' + "^B,H/)- the general case in which we do not follow 
the approximation given in eq. (2.3), the inelasticity, and with it the attenuation 
length, may be dramatically different from the mean free path, and therefore has to 
be considered in more detail. 

4.3 Spontaneous Decay 

Spontaneous decay of nuclei is a process described by -> ^',N' + ^N". In this 
work we are mainly interested into decays in which there is one heavy daughter 
nucleus and possibly other much lighter fragment (s). Some nuclei can decay spon- 
taneously even without LIV. This occurs mainly for nuclei with either proton or 
neutron excess emitting (most probably) single nucleons (e.g. 2He — )■ gHe + qU or 
gB — 4Be + \p) and is possible if the resulting nuclei have a total mass smaller than 
the mother nucleus one. 

In the presence of LIV, however, spontaneous decay may happen even for nuclei 
which cannot spontaneously decay in the LI framework. The main difference with 
respect to the LI case is that now there is a certain energy threshold for the decay 
to occur. To analyze the kinematics of the process we use eq. (4.4) by setting e = 0, 
(i.e. there is no low energy photon to interact with) since with this substitution the 
energy-momentum conservation eq. (4.3) still holds. The result is 



This equation has either none or exactly one real and positive solution for p depending 
on y. In the former case there is no spontaneous decay, similar to the LI situation, 
while in the latter case the (lower) spontaneous decay threshold p^^^ is obtained by 
minimizing the result with respect to y. 

Since the nuclei considered become unstable under such conditions, this kind of 
decay can be assumed to take place on a very short time scale. ^ This would mean 

■^In the absence of a proper theory of nuclei disruption it is impossible to reliably estimate this 
timescale. However, once above threshold this decay process is similar to the standard spontaneous 
decay, whose general timescales range from fractions of seconds to days or years. Moreover, the 
decay phase space grows with p'* if n = 2, hence the lifetime is expected to shorten with increasing 
energy. 




P' 



+2 = Mi 



(4.10) 
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that such particles cannot be observed on Earth. Analyzing the highest energy events 
under these aspects would therefore give information about the values of the LIV 
parameters. 

4.4 Vacuum Cherenkov emission 

The Vacuum Cherenkov effect (VC), i.e. the emission of a photon by a superlumi- 
nal charged particle in vacuum (in the case here ^V — )■ + 7) is, just like the 
spontaneous decay described above, forbidden for the Lorentz invariant case due to 
energy-momentum conservation. In the LIV framework, however, it becomes possible 
and must be taken into account [33, 20]. 

Since the reaction time scale for VC is usually very small [1], as low as Ins 
in the case of lOTeV electrons, the mean free path for the reaction is negligibly 
small compared to intergalactic propagation. Therefore, if a nucleus is observed 
with a certain energy the VC energy threshold must be larger than E. Moreover, 
assuming that LIV is much smaller for photons than for nuclei [22, 23, 24], we consider 
only emission of soft photons. This is a valid assumption, as the phase space of the 
reaction just above threshold is large enough to warrant the short mean free path 
necessary to set constraints [20]. 

We compute again the threshold momentum by exploiting energy-momentum 
conservation 



[E'^^z + eyc) -{Pa,z+P,,vc) . (4.11) 

After performing a calculation similar to the one done for eq. (4.3), the momentum 
threshold is [33] 

n+2 _ <,zM^l 



which leads to the constraint 



- ^ (A 12) 



if a nucleus with energy ~ E' is observed. 



5. Results 

We show here results obtained in the simplified case described in sec. 2.1, considering 
n = 2 with Vaz~a^^V being the LIV coefficient of a nucleon, and focusing on single 
nucleon emission. 

Assuming that current hints for a heavy composition at energies E ~ 10^^'^ eV 
[34] may be confirmed in the future, and that some UHECR is observed up to ~ 
10^° eV [52], one could place a first constraint on the absence of spontaneous decay 
for nuclei which cannot spontaneously decay without LIV as well. It will place a 
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Figure 2: Threshold momenta of spontaneous decay (left) and VC emission (right) for 
different elements. The respective shaded regions are excluded when the corresponding 
nucleus is observed at a given momentum value. Two of such observations are marked by 
the thick horizontal lines, taken at approx. lO^^'® eV from [53] and at 10^*^ eV from [52]. 

limit on ?7 < 0, because in this case the energy of the emitted nucleon is lowered 
with respect to the LI case until it "compensates" the binding energy of the nucleons 
in the initial nucleus in the energy- momentum conservation. We show in fig. 2 (left 
panel) the behavior of the threshold momentum with respect to rj < for three 
representative nuclear species: ^^Fe, ^^0 and ^He. If the UHECRs at the highest 
energies were Fe, then the limit is r/ > —1 for 10^^'^ eV (and —3 x 10~^ for 10^'' eV), 
while if they were lighter species, like He, then a remarkable constraint rj > —3 x 10~^ 
(—7 X 10~^) could be placed. An upper limit for rj > can instead be obtained from 
VC. This is illustrated in fig. 2 (right panel). Assuming UHECR to be mainly iron 
at the highest energies the constraint is given by r/ < 2 x 10^ for nuclei observed at 
lO^^-^ eV (and < 4 for lO^o eV), while for helium it is r/ < 4 x IQ-^ (10"^). 

A distinct upper constraint could in principle be obtained by assuming that pho- 
todisintegration is at work in UHECRs and then computing the largest value of r] 
for which the reaction is kinematically allowed in the required energy range. The 
change of the threshold momentum due to LIV is, however, very small as it can 
be understood with the following reasoning: The threshold equation (4.4) gives a 
situation similar to the LI case for small r), i.e. the threshold inelasticity is approxi- 
mately ythr ~ UthrLi = ' • For large positive r] the first term in eq. (4.4) 

mA',z' + rriB^w 
becomes the dominating one, so the equation reduces to 

r/ (1-2/)"+^^ ^ f^O. (5.1) 

Solving this gives Uthr- For single nucleon emission it is ythr ~ \ (which can be 
confirmed by plugging in), while for larger B at ythr ~ ^ the left hand expression of 
eq. (5.1) is at least as close to zero as possible. Substituting this back into eq. (4.4) 



and solving it for p gives the lower threshold for large t]. For 5 = 1, 

hmp^, ^ — ( —. — ^+"^iM/ • (5-2) 

We checked a posteriori that the approximation, eq. (5.1), is valid. Compared to the 
LI case P^^^j^j = ^ {{^A' + raif — m^), the relative difference is 

^Pthr ^ ^thr ~ hhr,LI _ ( ~ "^0 ^-^ 

Pthr,LI Pthr,LI + niif - m\ 

which is independent of 77 in the limit of large 77. For different elements we obtain 
= 7 X 10"^ for ^He, 2 x lO"-'^ for ^^O and 3 x 10"^ for ^^Fe. Hence an upper 

?-thr,LI 

constraint cannot be set just by considering the threshold shift because the relative 
change does not exceed 3 x 10"^ for elements up to iron and is therefore too small 
compared to the accuracy required to match the data. 

On the other hand, statistical quantities related to the probability of the reaction 
to happen, i.e. the different propagation lengths, can lead to interesting results. Now 
the question might arise how it is possible that these quantities are significantly 
affected by LIV while the thresholds are not. Indeed, the photon energy shift in the 
CR rest frame is linear with respect to 77, thus making the phase space and with it 
the cross section at a given photon energy e be dramatically affected by LIV. 

The mean free path for different elements and for different values of 77 is shown 
in fig. 3. LIV introduces quite interesting deviations from the LI case. The mean free 
path differs by at least one order of magnitude at the highest energies. The deviation 
goes however in different directions: For positive LIV parameters it decreases with 
growing momentum down to distances below 1 Mpc and then has a cutoff due to 
VC. On the other hand, for negative values it first increases (because of the reduced 
phase space), and then ends with a cutoff at the spontaneous decay threshold, i.e. a 
very fast decrease of Xmfp down to cosmologically tiny values. 

This produces interesting consequences on the UHECR spectrum. First of all, 
assuming that the closest UHECR source lies at least at a few Mpc distance from 
the Earth, the UHECR spectrum should display a sharp cutoff of the heavy nuclei 
component corresponding to energies for which Xmfp ^ 1 Mpc. On the other hand, 
if UHECR nuclei are found at such large energies, with such a small MFP, then a 
UHECR source should be present at a distance closer than 1 Mpc from the Earth. 
This would produce a strong anisotropy of the flux towards the source direction at 
E > 10^^'^ eV. As it is clear from fig. 3, this reasoning has the potential to set a 
constraint |7/| < 10~^ even for ^^Fe, once better data will be available. Given that the 
experimental data at such high energies still suffer from statistical and systematic 
uncertainties, we leave for future work the quantitative determination of the LIV 
UHECR spectrum and anisotropy, for which a full simulation of the propagation 
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Figure 3: The mean free path of photodisintegration in the single parameter model for 
s^Fe (top), (middle) and ^He (bottom). 



of heavy nuclei in the IGM, following the complete photodisintegration pattern, is 
needed. 
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E^,, = 1020eV 


^He 


-3 X 10-3 < ^ < 4 X 10-3 


-7 X 10-5 < 7] < 1 X 10-^ 


160 


-7 X 10-2 < < 1 


-2 X 10-3 < < 3 X 10-2 


56pg 


-1<^<200 


-3 X 10-2 < r/ < 4 



Table 2: Constraints given by spontaneous decay and vacuum Cherenkov emission if a 
nucleus of the species given in the left column with energy E^ax is observed. 

6. Conclusions 

Motivated by recent experimental results hinting at a significant presence of heavy 
nuclei in UHECRs, we analyzed the effects of a possible Lorentz invariance viola- 
tion on the propagation of ultra high energy cosmic ray nuclei in the intergalactic 
medium. Relevant effects were found in the process of photodisintegration onto in- 
tergalactic medium background photon fields. Moreover, we studied also Lorentz 
invariance violation effects on the nuclear spontaneous decay, and on possible emis- 
sion of Cherenkov radiation by superluminal charges in vacuum. While we did the 
computation of the main effects for a general modified dispersion relation, we ob- 
tained results only in the case of 0(£'2/Mpj) modified dispersion relations, assuming 
also that there is only one Lorentz invariance violation parameter controlling Lorentz 
invariance violation for all nuclear species, rj. 

A lower limit on r] can be set by requiring that nuclei, which are stable in 
the Lorentz invariant case, should not undergo spontaneous decay emitting single 
nucleons below energies of the order of lO^^'^ to 102*^ eV, at which a substantial 
fraction of heavy ultra high energy cosmic ray nuclei has been observed while an 
upper limit is given by an analogous argument for vacuum Cherenkov emission. 

Constraints obtained with these two techniques are summarized in Tab. 2. These 
constraints are to be compared to the results obtained assuming ultra high energy 
cosmic rays to be only protons [20]. In [20] it was found, for the case of pure 
proton composition, — 10^^ < 77 < 10"^. As it can be seen, our constraints improve 
when lighter nuclei are considered, and agree well with the ones presented in [20] 
when extrapolated to pure proton composition. We remark however that the recent 
findings of the presence of heavy nuclei in the ultra high energy cosmic ray spectrum 
might in fact invalidate the hypothesis of pure proton composition adopted in [20]. 

It is notable that we did not derive bounds coming directly from Lorentz invari- 
ance violation effects in photodisintegration. On the one hand, present experimental 
data on the chemical composition are still affected by systematics. On the other hand, 
to actually compute the theoretically expected spectrum and chemical composition 
in the presence of violations of Lorentz invariance requires strong computational ef- 
forts. Therefore, at present it is not possible to derive robust limits from the spectral 
information. Nevertheless, we derived expected order-of-magnitude limits assuming 
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that photodisintegration be at work at energies of order lO^^'® eV for all nuclear 
species, and we leave a more detailed evaluation of the constraints for future work. 
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